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Abstract 

In this paper, we investigate large amplitude solutions to a system of conservation laws 
which is transformed, by a change of variable, from the well-known Keller-Segel model de- 
scribing cell (bacteria) movement toward the concentration gradient of the chemical that is 
consumed by the cells. For the Cauchy problem and initial-boundary value problem, the 
global unique solvability is proved based on the energy method. In particular, our main 
purpose is to investigate the convergence rates as the diffusion parameter e goes to zero. It 
is shown that the convergence rates in L°°-norm are of the order O (e) and Oie^l^^ corre- 
sponding to the Cauchy problem and the initial-boundary value problem respectively. 
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1 Introduction 

In this paper, we investigate large amplitude solutions to a system of conservation laws 



X 



(1.1) 



which is transformed, by a change of variable, from the well-known Keller-Segel model describing 
cell (bacteria) movement toward the concentration gradient of the chemical that is consumed 
by the cells. Here £ > is a positive constant. 

The aim here is to study the global unique solvability on the Cauchy problem and initial- 
boundary value problem in the framework of large-amplitude solutions. In particular, our 
main purpose is to investigate the convergence rates as the diffusion parameter e goes to zero. It 
is shown that the convergence rates in L°°-norm are of the order O (e) and 0{e^^^) corresponding 
to the Cauchy problem and the initial-boundary value problem respectively. 

Firstly we are concerned with the Cauchy problem of (1.1) with initial data 

{u^,v^) (x,0) = {uI{x),Vq{x)) {Q,Voo) as X ^ ±oo. (1.2) 

Hereafter, v^o is a given positive constant. On one hand, the global existence of solutions to the 
Cauchy problem (1.1), (1-2) for any fixed e > is shown. On the other hand, we prove that 
solutions to the Cauchy problem (1.1), (1.2) convergence to solutions to the problem of zero 
viscosity limit when e — > 0+, which is formally formulated as follows: 



u^t-Vx = 0, 



(1.3) 

X ~ "XXI 



with initial data 

{x, 0) = (i4{x),vl{x)) ^ (0, v^) as x^ ±oo. (1.4) 
In addition, we consider the initial-boundary value problem of (1.1) with initial data 

{v'y){x,i)) = {ul{x),vl{x)), 0<x<l, (1.5) 
and the boundary conditions 

(n^<)(0,^) = (n^<)(l,^) = (0,0), i>0, (1.6) 

which implies 

<,(0,t) =nL(l,0 = 0, t>0, (1.7) 

where the compatibility conditions tio(O) = ^io(l) = and Vq^[{)) = Vq^{1) = arc satisfied. 

For one thing, the global existence of solutions to the initial-boundary value problem of 
(1.1), (1.5), (1.6) for any fixed e > is shown. For another thing, we prove that solutions to the 
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initial-boundary value problem of (1.1), (1-5), (1.6) convergence to solutions to the problem of 
zero viscosity limit when e ^ 0'^, which is formally formulated as follows: 



t ~ - 0' 



with initial data 



{u'^y) (x,0) = {ul{x),vl{x)) , < X < 1, 
and the boundary conditions 

tt°(0,i) = it°(l,t) = 0, i>0, 



(1.8) 

(1.9) 

(1.10) 
(1.11) 



which implies 

i;°(0,i) = ^;0(l,i) = 0, t>0, 

where the compatibility conditions Mo(0) = Wo(l) = and Vq^{^) = ^^^^^(l) = are satisfied. 

Precisely speaking, these results are given as follows. We have the first main theorem, which 
is concerned with the global existence of the Cauchy problem (1.1)-(1.2). 

Theorem 1.1 (Global existence) Let {uq,Vq — Voo) G H'^(R) and assume thai there exists a 
positive constant a > such that in^t;Q(a;) > a. Then there exists a unique global solution 
{u^{x,t),v^{x,t)) of the Cauchy problem (1.1)-(1.2) which satisfies 

' (n^ V' - v^) e L°°([0, oo), iJ2(R)), 

^ e L2([0,oo),L2(R)), v^^, G L\[0,oo),H\K)) 



(1.12) 



and 



rt 

\u'{t)\\H2 + \\v'{t)-Voo\\H2+S / \\ul{s)fH2ds 

Jo 



+ 







+ 



L2 



+ \\vU(s)\\'ii > ds 



(1.13) 



< 



' i\-'XX\"J\\H'^ ^ 
\Uq\\h2 + 11^0 - '"oo\\H2^ 

Moreover, for any fixed T > 0, there exists a positive constant C (T) which depends only on T 
such that 



v%x,t) > C(T). 
Here C and C{T) are all positive constants independent of e. 



(1.14) 



Remark 1.1. In the case of e = 0, the global unique solvability on the Cauchy problem (1.3)- 
(1.4) was included in [5], which will be summarized in Lemma 2.4 later. As is well-known that it 
is more difficult when the nonlinearity is higher in the setting of large amplitude. It is nontrivial 
to obtain the estimate on (1.13) in the setting of large amplitude for the case of e > since 
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the nonhnear term appears. We overcome these difficulties by applying some suitable 

interpolation inequalities. 

Next, the following main result refers to convergence rate under the additional assumptions 
on initial data. 

Theorem 1.2 (L°°-convergence rate) Under the same conditions of Theorem 1.1, we as- 
sume the initial data {uq,Vq) and (^Uq,Vq) satisfy \\uq ^ 'f^olli/2(R) 
the further regularity [uq,Vq) G H^(R). Then we have 



+ H 



V, 



01 



(u' - (t) 



+ 



L°°(R) 

Here C > is a positive constant independent of e. 



(v' - v'^) (t) 



L°°(R) 



< Ce. (1.15) 



Remark 1.2. In the case of e > 0, for the initial-boundary value problem of (1.1) with Dirichlet 
boundary conditions u^{{),t) = u^{l,t) = 0, v^{0,t) = v^{l,t) = 0, with the help of the methods 
in [25], we can prove the global unique solvability in the setting of small amplitude since the 
difficulties derived from the boundary effect on the viscosity eul^. What's more, the boundary 
layer will occur and the boundary layer thickness will be shown as in [4, 7, 18]. These are left 
to our future study since these are not in accordance with our main aim in the paper, which 
mainly refers to large amplitude solution. 

Our final result focusing on the initial-boundary value problem is stated as follows. 

Theorem 1.3 (Initial-boundary value problem) Let {uq,Vq) G i?^([0, 1]) and assume that 

there exists a positive constant a > such that inf Vn(x) > a. 

xe[o,i] ~ 

(i) Then there exists a unique global solution (u^(x,t),v^ {x,t)) of the initial-boundary value 
problem (1.1), (1-5) and (1.6) which satisfies 

^ ^ G L2([0,oo),L2([0,1])), G L2([0,oo),i7H[0,l])) 



(1.16) 



and 



\u%ml, + \\v^{t)\\%.+e f\\uUs)fH.ds 

Jo 



+ 



f 

Jo 



+ 



L2 



V^<{s)\\^, + \Kx{s)\\m\ds 



(1.17) 



|2 , II e||2 



<c(ll^ol 

Moreover, for any fixed T > 0, there exists a positive constant C{T) which depends only on T 
such that 



v^{x,t) > C{T). 



(1.18) 



(ii) We assume the initial data (uq, Vq) and (uq, Vq) satisfy ||no — ^ollij2([o i])"'"ll'^o ~ '^ollij2([o i]) 
0{e). Then we have 



{u^ - u') (t) 



+ 



(v' - ?;°) (t) 



L°°([0,1]) 



(1.19) 
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Here C and C{T) are all positive constants independent of e. 

Remark 1.3. Convergence rates (1.19) on the initial-boundary value problem is slower than 
one in (1.15) on the Cauchy problem, which is caused by the boundary effect. Essentially it 
is different between initial-boundary value problem and Cauchy problem that the regularity of 
solutions to the Cauchy problem can be improved when the regularity of initial data is imposed 
further regularity whereas it fails for the initial-boundary value problem. 

The problem of the zero viscosity limit is one of the important topics. In particular, when 
parabolic equations with small viscosity are applied as perturbations, convergence rate of the 
Cauchy problem or the boundary layer question of the initial-boundary value problem for many 
other equations also arises in the theory of hyperbolic systems in the case of one-dimension or 
multi-dimension, cf. [3, 4, 7, 18, 22, 23]. To our knowledge, fewer results on the equations (1.1) 
have been obtained in this direction. 

Now let us review some known results related to the system (1.1) and (1.3), which has been 
extensively studied by several authors in different contexts, cf. [5, 11, 14, 25]. The conservation 
laws (1.1) are derived from the original well-known Keller-Segel model 

Ct = eCxx - Uf{c), 

(1.20) 

which was proposed by Keller and Segel in [8] to describe the traveling band behavior of bacteria 
due to the chemotactic response (i.e., the oriented movement of cells to the chemical concentra- 
tion gradient) observed in experiments [1, 2]. In model (1.20), u{x,t) and c{x,t) denote the cell 
density and the chemical concentration, respectively. D > is the diffusion rate of cells (bacte- 
ria) and £ > is the diffusion rate of chemical substance, x is a positive constant often referred 
to as chemoscnsitivity. /(c) is a kinetic function describing the chemical reaction between cells 
and the chemical. 

When /(c) is a positive constant, namely, /(c) = a > 0, the existence of traveling wave 
solutions of (1.20) with e = was established by Keller and Segel themselves in [8]. When 
e 7^ 0, the existence and linear instability of traveling wave solutions of (1.20) were shown by 
Nagai and Ikeda in [15] where the authors also obtained the diffusion limits of traveling wave 
solutions of (1.20) as e approaches zero. Precisely they proved that the traveling wave solution 
in the form {B^{x,t),S^{x,t)) = {B^ (z) , {z)) {z = x — ct) approximates to the corresponding 
to traveling wave solution (B'^{x,t), S^{x,t)) = (B^ {z) , {z)) {z = x — ct) when e goes to zero. 
For the reduction of Keller-Segel system (1.20) to system (1.1), we refer to Refs. [14], see also 
Appendix in this paper. 

In [14], Li and Wang established the existence and the nonlinear stability of traveling wave 
solutions to a system of conservation laws (1.1). They prove the existence of traveling fronts by 
the phase plane analysis and show the asymptotic nonlinear stability of traveling wave solutions 
without the smallness assumption on the wave strengths by the method of energy estimates. 

There is limiting case of the Keller-Segel model (1.20). That is when the diffusion of chemical 
substance is so small that it is negligible, i.e, e 0+, then the model (1.20) becomes 
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f 

ct = -uf{c), 

(1.21) 

A version of system (1.21) was proposed by Othmer and Stevens in [16] to describe the 
chemotactic movement of particles where the chemicals are non-diffusible. Othmer and Stevens 
in [16] have developed a number of mathematical models of chemotaxis to illustrate aggregation 
leading (numerically) to nonconstant steady-states, blow-up resulting in the formation of singu- 
larities and collapse or the formation of a spatially uniform steady state. The models developed 
in [16] have been studied in depth by Lcvine and Slccman in [9]. They gave some heuristic un- 
derstanding of some of these phenomena and investigated the properties of solutions of a system 
of chemotaxis equation arising in the theory of reinforced random walks. Y. Yang, H. Chen and 
W.A. Liu in [24] studied the global existence and blow-up in a finite-time of solutions for the 
case considered in [9] , respectively. They found that even at the same growth rate the behavior 
of the biological systems can be very different just because they started their action in different 
conditions. For the other results on the initial-boundary value problem of (1.21) refer to [6, 10]. 

Recently, the modified model related closely to chemotaxis was also investigated by Painter 
and Hillen in [17]. They explored the dynamics of a one-dimensional Keller-Segel type model 
for chemotaxis incorporating a logistic cell growth term. 

Finally, we have to mention the work in [25, 5], which investigated the initial-boundary 
value problem and Cauchy problem on limit equation (1.3) and motivate our investigation in 
this paper. Zhang and Zhu in [25] studied the initial-boundary value problem of (1.3) with 
initial data 

(^o,t;0)(x,0) = (u[l(x),^;0(x)) (1.22) 

and boundary conditions 

u'^{0,t) = u^{l,t) =0. (1.23) 

Let [uq,Vq — 1) G //^([O, 1]) and assume that there exists a positive constant e > suffi- 
ciently small such that ||^o||/f2 + H'^q — l||^2 < £• Then there exists a unique global solution 
(vP{x,t),v^{x,t)) of the initial-boundary value problem (1.3), (1.22) and (1.23) which satisfies 

(u'^y) G L-([0,oo),//2([o,i])), g L-'{[0,oo),H\[0,l])). (1.24) 
Guo etc. in [5] studied the Cauchy problem of (1.3) with initial data 

(x,0) = (u'^{x),v^{x)'^ iu±,Voo) as x ^ ±oo. (1.25) 

They proved the existence of global solutions to the Cauchy problem of a hyperbolic-parabolic 
coupled system with large initial data, which is summarized in Lemma 2.4 later. 

Notations: Throughout this paper, we denote positive constants by C may depends on T, but 
is independent of e. Moreover, the character "C" may differ in different places. LP = LP{Q) 
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(1 < P < oo) denotes usual Lebesgue space with the norm 



11/11 



LP(Q) 



1 < p < oo 




sup|/(a;)|. 



(Q) (/ > 0) denotes the usual /th-order Sobolev space with the norm 



( 



1 

2 



11/11 



E 11^^/11' 



where i7 = R or [0,1], and || • || = || • ||o = || • Hl^- For simplicity, ||/(-,t)||LP and ||/(-,i)||i are 
denoted by ||/(t)||LP and ||/(t)||z respectively. 

The rest of the paper is organized as follows. In Section 2, wc study the Cauchy problem 
(1.1) and (1.2). For any fixed e > 0, the global unique solvability are obtained. When diffusion 
parameter e — t- 0"*", the limit problem is considered and we show that the convergence rates in 
L°°-norm is of the order 0(e). In Section 3, we study the initial-boundary value problem (1.1), 
(1.5) and (1.6). In the last Section 4, for the convenience of the readers, we use the appendix 
to give the derivation of system (1.1) and (1.3). 

2 Cauchy problem 

2.1 Global existence on the case of e > 

In this section, wc arc concerned with the global existence of large-amplitude solutions to 
the Cauchy problem (1.1), (1-2), which will be proven by continuing a unique local solution 
using a uniform-in-e a priori estimate in the setting of large data. The construction on the local 
existence of the solutions is standard based on iteration argument and Fixed Point Theorem. 
Next, we will devoted ourself to obtaining some uniform-in-e a priori estimates on the solution 
{u^{x,t),v^{x,t)) to the Cauchy problem (1.1), (1.2) for all £ > 0. To this end, v^{x,t) is 
supposed to satisfy the priori assumption 



for any fixed T > 0. 

Lemma 2.1 (Basic energy estimate) Assume that the assumptions listed in Theorem 1.1 
are satisfied. Then there exists a positive constant C , independent of e, such that 



v^{x,t)>{), (x, t) G R X [0, T] 



(2.1) 




dxds + £ 




I / ('"D^ dxds = / rjodx 
Jr 



(2.2) 
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Proof. As in [25, 5], it is easy to see that 



q v^) = -u'v' In — + 1^ {u'f 



is an entropy-entropy flux pair to the hyperbohc system (1.1) normahzed at (0,Woo)- Moreover, 
based on the definition of the entropy-entropy flux pair (see [20]), it is easy to verify that 
(77 (u^, v^) ,q{u^ ,v^)) satisfies 



rit + (lx = + v^x In — • 



(2.3) 



Integrating the above identity with respect x and t over R x [0, t], by using integration by parts 
we can get (2.2). This completes the proof of Lemma 2.1. 

Lemma 2.2 (First order energy estimate) Assume that the assumptions listed in Theorem 
1.1 are satisfied. Then there exists a positive constant C , independent of e, such that 



lL°°(Rx[o,r]) 



+ 



Il-(Rx[o,t]) + / {{ulf + {ylf\ dx 

{v' {u%f + {v%,f} dxds + e f j {u%,f dxds < a 



rt 

Here and hereafter, C denotes positive constant depends on T. 



(2.4) 



Proof. As in [5], notice that 



It. 



xt '^xx 



(2.5) 



Multiplying (2.5) by 2ii§, then integrating the resulting equation with respect x and t over 
R X [0, t], we have by exploiting some integrations by parts and Cauchy-Schwarz inequality that 

/ {ulfdx-j {ul^fdx + 2f f v'{u%fdxds + 2e f j {u%^f dxds 
Jn Jr Jo Jr Jo Jr 

= 2 / / v^u^dxds — 2 / / u^v%u%dxds — 2e / (u|)^ dxds 
Jo Jr Jo Jr Jo Jr 

= j2h. (2.6) 
1=1 

/1-/3 are estimated as follows: 

h = —2 [ {vq — Voo) UQ^dx + 2 [ {v^ — Voo) u%dx — 2 / f {v^ — Voo) u^^dxds 
Jr Jr Jo Jr 

= -2 [ (wq - Voo) UQ^dx + 2 [ (v^^ - Voo) u%dx - 2 [ / (f ^ - Voo) v%^dxds 
Jr Jr Jo Jr 

+2e / (v^ — Voo) {U'^)'^ dxds — 2e / {v^ — Voo) u^r^r^dxds 

Jo Jr I ixx Jq 

= E^i- (2.7) 

i=l 
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It is easy to get 

II < [ {v'o- v^f dx+ I {ulJ" dx < C. (2.8) 
Notice that there exists a positive constant C > such that 

-Voo\< civ' hi— -{v^ -Voo)\ <Cr]{u',v'), iovv'eQ.i (2.9) 



and 



since 



and 



-Voo\^ <c[v'hi— -{v' -Voo)\ <Cri{u',v') , for G O2, (2.10) 



- 'fool 



v'hi^ - {v^ - Voo) 3 In I - 1 



^6^2 v^ln^-{v' -Voo) ~ 61n|-2' 



Here 
and 



^1 = {v : \v — Uoo < t = {v : — — <v < 



2 J I 2 - - 2 
Then we have from (2.9) and (2.10) 

/i < 2 / {v' - v^f dx + l [ {ulf dx 

= 2/ {v' -Voofdx + 2 [ (v' - Voof dx + I [ {u%fdx 

< C (1 + \\v' - ?^oo||l-(Rx[o,t])) mdx + \J^ (Kf dx. (2.11) 
By using integration by part, Cauchy-Schwartz inequality and Lemma 2.1, we have 

If = 2 /* / Kfdxds < ||i;^IL^(Rx[o,T]) r / ^dxds, (2.12) 

Jo t/R Jo t/R " 

if = —As [ [ u'u^v^dxds 
Jo Jr 

< ll«1lioo(Rx[o,T]) / / -^dxds 
Jo «/R, Jo »/R " 

< f f vH<fdxds + C\\u^l^^^^^^^j^^ (2.13) 

Jo J R 
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and 

if = 2e [ [ vl.u%^dxds 
Jo JR 

I {<x? dxds + e\\v'\\Loo^^y,yQ^T]) I I '^d^ds. (2.14) 

I' R. ^ ^ R ^ 

The estimate on I2 is the same as those of if. 

Finahy, we estimate I^. By using GagUardo-Nirenberg inequahty and Young inequahty, it 
fohows that 

||<(t)||i3 <C||<,(t)||i U{t)\\i <\\\u%,{t)\\^ + C\\u%t)r. (2.15) 
Thus, we have from Lemma 2.1 

h<^£j^J {ul^f dxds + Ce j^^i^j {u^fdx^ ds 

<\ef ! {ul,fdxds + C. (2.16) 
Substituting (2.7)-(2.16) into (2.6), we get that 

(u^.)^dx+ / / f ^ (tt|)^ dxds + e / / {u^.^)'^ dxds 



R Jo JR Jo JR 







1 Kfdx^ 




JR 


Jo Jr 



< C{T) (1 + ||'u^||ioo(Rx[o,T]) + II'"^IIl°°(Rx[o,t])) ■ (2-17) 
Prom the Sobolev inequahty and Cauchy inequahty , we have 

ll«'llioc(R,[o,T]) < 2 sup {||u^(t)|| \\ul{t)\\} < sup \C \\u%t)f + ^ \K{t)f] . (2.18) 

Substituting (2.18) into (2.17), using Lemma 2.1, we deduce immediately 

it|)^ dxds + e / {u^x)^ dxds 
Jo Jr 

<c(l + lb'llLoc(Rx[o,r]))- (2-19) 
Prom Lemma 2.1, (2.18) and (2.19), we have 

II'"^IIl°°(Rx[o,t]) + ||^'^||x,<x.(rx[o,t])) • (2.20) 

The rest of the proof of this lemma is similar to those in [5]. By differentiating the second 
equation of (1.1) with respect to x once and multiplying the resulting identity by 2v^, then 
integrating the final equation with respect x and t over R x [0, i], we have by integrations by 
parts, (2.18), (2.19) and Lemma 2.1 that 

{vlf dx + fj^ {v%,f dxds <c(l + ||^^||ioc(Rx[o,r])) • (2-21) 
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In order to get the L°°(R x [O,r])-norm estimate on {u'^{x,t),v^{x,t)), take 

e \ I Jh{z)dz, for >Voo, 
h{v') = v'hi {v'-v^), ^{v')=l-'-^^ (2.22) 

[ j yjh{z)dz, for < < v^. 
Prom (2.9) and (2.10), there exists a positive constant C > such that 

\V' -Voo\<C{h{v') + l), 

which implies 



l^^'-t-ool^ <C[Jh{v^) + l 



(2.23) 



Prom (2.22) and (2.23), we have 



and 



\ {v^ -Voo)^ = r \z-Voo\^dz<C r Jh{z)dz + C{v^ - v^) , for > (2.24) 

-^{voo-v^)^= \z-Voo\^dz<C J h{z)dz + C {voo - v^) , for < i;oo-(2.25) 
On the other hand, we have 



ip{v%dy= / Jh{v^Xydy 

-OO J —CO 

\\vl\\ <c(l + \\v^\\l'>°{-rx[o,t])) ■ 



Thus for any > and some positive constant C > 0, we have 



(2.26) 



c(l 



+ ll^^llL°°(Rx[0,T])j • (2.27) 



This means that 



I^^IIl°°(Rx[0,T]) ^ ^' 



and consequently we have from (2.20) 



l'"^llL°°(Rx[0,r]) — ^■ 



(2.28) 



(2.29) 



Combination of (2.19), (2.21), (2.28) and (2.29) yields (2.4). This completes the proof of 
Lemma 2.2. 

Remark 2.1. The assumption {uq,Vq — Vqo) € i?^(R) implies 

j^riodx < C (ll^xgf + \\vl - v^f) 
which and (2.4), (2.11) show that 

\\v'{t) - v^f < C {\\ul\\^ + \\vl - v^\f) . 
At last, we deduce the L°°([0,T],L^(R))-norm estimate on {ul.^{x,t),vl.^{x,t)). 
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Lemma 2.3 (Second-order energy estimate) Assume that the assumptions listed in Theo- 
rem 1.1 are satisfied. Then there exists a positive constant C , independent of e, such that 

+ dx + {e {u%,,f + {v%,,f] dxds < C. (2.30) 

Proof. Differentiating (1.1) with respect to ,t, multiplying the resulting identity by — 2u|^2, and 
—2v%,j.,j. respectively, and integrating the adding result with respect x and t over R x [0,t], we 
have from some integrations by parts that 

+ dx + 2 {e + dxds 

= {{uoxxf + (voxxf) ~ 2 ^ dxds 

+^ f f ^'xxxKxdxds-2 f I {e{u^f] u%Jxds 
7 

= E^- (2.31) 

i=4 

The estimates on I^-Iq can be found in [5] as follow: 

h + h< f f ivL,f dxds + c\l+ f [ ^^^A (2.32) 

Jo Jr I Jo Jr ) 

Next we will be devoted to estimating iy. 

I7 = 4£ [ [ {ulf ul^^dxds + Ae ( f ul^u%^^dxds 
Jo Jr Jo Jr 

= Y,I\. (2.33) 

1=1 

From Cauchy inequality and Lemma 2.1, one has 



' I I {'^Ixxf dxds + 16e f f «)^ dxds 
Jo Jr Jo Jr 

< \e J' dxds + Cj^y^ [u^^ dx^'ds 

<\^f I {u%,,fdxds + C. (2.34) 
2. Jo Jr 

Here we used the fact obtained by Gagliardo-Nirenberg inequality and Young inequality 

||<(t)||i4 <ChL,(t)||i \\u^{t)\\l < l||nL.(i)f + C||u^(t)f^ (2.35) 

In addition, one has by Cauchy inequality and Lemma 2.2 

l7<\^ J^iKxxf dxds + 16s Wu^'Wl^^-^^^o^T^y^ j^{u%^fdxds 

<\e f I {ul,,fdxds + C f [ {u%,fdxds. (2.36) 
4 Jo Jr Jo Jr 
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Substituting (2.32)-(2.36) into (2.31), we have 

<c[l + fj^ dxds^ . (2.37) 

(2.30) immediately follows by applying the Gronwall inequality. This completes the proof of 
Lemma 2.3. 

Remark 2.2. By Sobolev inequality, Lemmas 2.2 and 2.3, we have 



KIIl-(Rx[o,t]) < Csup {|K(t)||^ II^LWIlH < c 



and 



II<IIl-(Rx[o,t]) <Csup{|K(t)p |bL(i)||H <C, 

which will be used later. 

Combination of Lemmas 2.1-2.3 and Remark 2.1 yields to the uniform-in-e a priori estimate 
(1.13), which implies the global existence of solutions by combining the local existence and 
uniqueness. Similar to [5], by standard maximal principle, we can prove (1.14) and the priori 
assumption (2.1) holds. This completes the proof of global existence in Theorem 1.1. 

2.2 Convergence rate of zero viscosity limit 

In this subsection, wc turn to our another result, which is concerned with convergence rates of 
the vanishing diffusion viscosity. That is, we will give the proof of Theorem 1.2, and it suffices 
to show the following Lemma 2.6. 

First, the global existence of the H'^ solutions to the Cauchy problem (1.3), (1.4) was included 
in [5], which is summarized as follows. 

Lemma 2.4 (Global existence on the limit problem, cf. [5]) Assume that the initial data 
{uq{x),Vq{x)) satisfies 

inf Vq {x) > a > and ( Uq — u, — ) G if2(R). (2.38) 

Then there exists a unique global solution {u^{x,t),v^{x,t)) of the Cauchy problem (1.3), (1-4) 
satisfying 

(u'^ -u,v° -vjj eL'^{[0,oo),H^(R)), w° G L2([0,oo),i^2(R)). (2.39) 

Moreover, for any fixed T > 0, there exists a positive constant C{T) > which depends only on 
T and \\ (uq — u,Vq — Voo) \\jj2 such that 

v^{x,t)>C{T). (2.40) 

Here the smooth monotone function u{x) satisfies u{x) = u±, ±x > 1 for fixed constants U- 
and u+ . 
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Next, we need improve the regularity on the solutions to the zero viscosity limit problem 
(1.3), (1.4). 

Lemma 2.5 (Regularity improved) Assume that the assumptions listed in Theorem 1.2 are 
satisfied. Then there exists a positive constant C such that 

[ + (^L.)'} dx + ^ {vl,,,)' dxds < C. (2.41) 

Proof. Differentiating (1.3) with respect to x three times, multiplying the resulting identity 
by li^xxx 2t;2^^ respectively, and integrating the adding result with respect x and t over 
R X [0, i], we have from some integrations by parts that 



R 



I (^Oxxx) ^ + (voxxx) ^ } ~ 2 / U^vlxx^lxxxdxds - ^ '^l^^x'^xxxxdxds 



— 6 / / u^^v^v^^^^dxds — 2 / / u^^^v'^v^^^^dxds + 2 
10 Jr Jo Jr 



t f ul 
Jr 



1=1 



J2-J& are estimated as follows: 

J2 <\ 1^(4^^^)" dxds + C 



2 



(2.42) 



%xxx ) dxds + C 



^ r / (vLY dxds, (2.44) 



J4<^J^J^ (v^xxx) dxds + C 



L°°(Rx[0,T]) 



J^(ul,f dxds, (2.45) 



2 



J5 + J6<IJ^ (vl^^^^dxds + C[l- 

Substituting (2.43)-(2.46) into (2.42), we get from Lemma 2.4 

{ ("Lj^ + (^Lx)^} {'^xxxxYdxds 
- ^ ^ /o7r (^°-)' ^^^^4 • (2-47) 



(2.41) immediately follows by applying the Gronwall inequality. This completes the proof of 
Lemma 2.5. 

Based on Theorem 1.1 and Lemmas 2.4-2.5, the following L^-convergence rates can be proved. 
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Lemma 2.6 (L^-Convergence rates) Assum,e that the assumptions listed in Theorem 1.2 are 
satisfied. Then there exists a positive constant C , independent of e, such that 

v2" 



dxdr < Ce 



(2.48) 



and 



R 



dx 



+ 



Jn 



dxdr < Ce 



Proof. Set 



V — V . 



(2.49) 



(2.50) 



Then we deduce from (1.1)-(1.2) and (1.3)-(1.4) that {x,t) satisfy the following Cauchy 

problem: 

>f + (£(n-)2-r)^ = eVL + e^L, 



with initial data 



(V^r)(x,o) = (V'S,0g). 



(2.51) 



(2.52) 



Proof of (2.48). 

Multiplying the fist and second equation of (2.51) by 2ip^ and 29^ respectively, integrating 
the adding result with respect x and t over R x [0, t] , we have 

^ + {e^f} dx + 2 1^ {s {^If + {eif} dxds 

= I {{i^lf + {dlf] dx - r / i)'u^u%dxds + 2e /* / i^^u^Jxds 
Jr JO JR Jo Jr 



-2 f I ^^e^dxds - 2 /* / i)^v^eidxds - 2 /* / u^O^Oldxds 
Jo Jr Jo Jr Jo Jr 
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i=7 

Here J7-J12 are estimated as follows. 
First, one has 



J7 < Ce^ 



(2.53) 



(2.54) 



since \\uf) - Mo||^2(r) + \\vf) - ^^o||h2(r) = 0(e). 

Next, from Cauchy-Schwarz inequality. Theorem 1.1 and Lemma 2.4 with u{x) = 0, we 
obtain 

■^8 < IK^|lioo(R^x[o^T]) / / {Ux)^dxds+ / / ((ip^)^ dxds 



{ip^f dxds, 



JR 



(2.55) 
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Jo Jn^ ' Jo Ju 

/* / {ip^fdxds (2.56) 
Jo Jr 



<Ce'^ + 



and 



+1 r / {eifdxds 

^ Jo Jr 

<C t i {{'^^f + {e'f] dxds + - I' I {e%f dxds. (2.57) 
Jo Jr^ ^ Jo Jn 

Substituting (2.54)-(2.57) into (2.53), we get 

1^ + (01'} dx + {e {i:lf + {eif] dxds 

< Ce^ + C j^{{'^^f + {0^?} dxds. (2.58) 

We immediately get (2.48) from Gronwall inequality. 
Proof of (2.49) . 

Multiplying the fist and second equation of (2.51) by — 2^0^^ and —29^^ respectively, inte- 
grating the adding result with respect x and t over R x [0,t], we have 

^ {m' + mf} dx +'^[1^ + dxds 

= 1 {ii'oxf + {0oxf}dx + 4e f I ijl^u'u%dxds-2e I f i^l^^^dxds 
Jn^ ^ Jo Jn Jo Jn 

-2/*/ r,J^dxds-2f [ ei^Jxds 
Jo Jn Jo Jn^ ' ^ 

17 
i=13 

Here J13- J17 are estimated as follows. 

First, similar to the estimate (2.54), one has 

^13 < Ce^ (2.60) 

Next, from Cauchy-Schwarz inequality. Theorem 1.1, Lemmas 2.4-2.5, Remark 2.2 and (2.48), 
we obtain 



Ji4 = -4e / / (u^.)^ dxds - 4e [ [ 
Jo Jn. Jo Jr 



ilil.u^ u%^dxds 

Rx[0,T]) 

{u%) dxds + Ce ||^i^||Loo(Rx[o,T]) {u%x) dxds 

+ [ [ {il^lf dxds 
Jo Jb, 



2 1 1 £ 1 1 2 



< Ce^ + /* / (O^ dxds, (2.61) 
Jo Jr 
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Ji5 -^^ {^Ix^ dxds 

Jo Jn 



+ / / [il^lfdxds 

JR 



dxds. 



(2.62) 



and 



Here 



Ji6 = 2 / / iljiei^dxds 

Jo JR 



< 



JR 



{rxfdxds + ] f f {ei^fdxds 

4 Jo JR 



Ji7 = -2 / / tp^v^ei^dxds -2 [ [ il^^viei^dxds -2 [ / ^t°^ 
Jo Jr Jo Jr Jo Jr 

-2 t I u^etei^dxds 



JR 



4 

1=1 



% 

17- 



1 /•* 



4 Jo Jr 



{OD dxds + C Wv^^W^oo^^^^^Q^T]) / / (V-x) c^arcfs 



JR 



<l/ / (eLfdxds + cf [ ii^lfdxds, 
4 Jo Jr Jo Jr 



(2.63) 



(2.64) 



(2.65) 



/.7a 

i^lxf dxds + C |K||ioo(Rx[o,T]) i'^^f dxds 



1 /•* 
< - 



>/R 



(2.66) 



JR 



L°°(Rx[0,T]) Jo Jr 



(r )2 dxds 



+C 



L°°(Rx[0,T]) 



/* / my 

Jo Jr 



dxds 



1 rt 
< - 



4 Jo Jr 



dxds + Ce . 



(2.67) 



Substituting (2.60)-(2.67) into (2.59), we get 

[irx? + {eif} dx + {e {^xxf + {0%,?] dxds 

KCs' + C 1^ /^{{^xf + mf} dxds. (2.68) 
We immediately get (2.49) from Gronwall inequality. This completes the proof of Lemma 2.6. 
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Finally, based on Lemma 2.6, we can prove Theorem 1.2. In fact, using Sobolev inequality, 
we also have from (2.48) and (2.49) 



< Ce 



(2.69) 



and 



L°°(R) 



< c 



(v^ - (i) 

< Ce. 

This completes the proof of Theorem 1.2. 



(2.70) 



3 Initial-boundary value problem 

In this section, we are concerned with the global existence of large-amplitude solutions 
to the initial-boundary value problem (1.1), (1.5) and (1.6). In particular, we investigate the 

convergence rates as the diffusion parameter e goes to zero. 

By the method similar to those in Subsection 2.1, we can prove (1.16) and (1.17) in Theorem 
1.3 (i). Now we devoted ourselves to claiming (1.18) in Theorem 1.3 (i). 
Proof of (1.18). From Sobolev inequality and (1.17), one has 



I'"IIIl°°([0,1]x[0,T]) ^ ^• 



Now set 



M= sup ||n^(x,t)||, t(;^(x,t) = e '^V(x,t), 
[o,i]x[o,T] 



where A > M is a positive constant. 

We can deduce from (1.1), (1.5) and (1.6) that w^{x,t) satisfies the following initial-boundary 
value problem 



wl-u'wl, + {X-ul)w' = wl,^, xeiO,l), Q<t<T, 
w^{x,0) = vf,{x) > a >0, xG[0, 1], 
wUO,t)=wUl,t)=0, 0<t<T. 
The standard maximal principle tells us that 

w%x,t) > 0, {x,t) G [0,1] X [o,r]. 



(3.1) 



(3.2) 



In fact, if w^{x,t) attains its negative minimum at some point {xo,to) G [0,1] x (0,r]. Thus 
w^{x,t) satisfies 



wl{xo,to) < 0, wl{xo,to) = 0, w^{xQ,to) < 0, wl^{xo,to) > 0, 
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which contradict to the first equation of (3.1) since X — > 0. 
One can deduce that w'^{x,t) satisfies from (3.1) and (3.2) 

wfix, t) + iX + M)w%x, t) - u'{x, t)wlix, t) > wl^ix, t). (3.3) 

Let H%x,t) = e-"* (e'^^+^^^w%x,t) - a), we have from (3.3) that 

HI{x,t) + {X + M)H%x,t)-u%x,t)H^{x,t)>H^^{x,t), xG(0, 1), 0<t<T, 

H^{x,0) =w^{x,0) -a>0, xe[0,l], 

i^|(0,i) = i?|(l,t)=0, 0<t<T. 
By the standard maximal principle, we have 

H%x,t)>0, (x,i) e [0,1] X [0,T]. 
Thus for any {x,t) G [0,1] x [0,T], v%x,t) satisfies 

v'{x,t)>ae~^'>ae-^^. (3.4) 

This completes the proof of (1.18). 

Similar results to those of Theorem 1.3 (i) can be obtained on the initial-boundary value 
problem on (1.8)-(1.11) as follows. 

Lemma 3.1 Let {uq,Vq) G i?^([0, 1]) and assume that there exists a positive constant a > 

such that inf Vq{x) > a. Then there exists a unique global solution {u^{x,t),v^{x,t)) of the 

xe[o,i] 

initial-boundary value problem (1.8)-(1.11) satisfying 
{u\v^)eL--{[Q,<^),H\[Q,l])), 

(3-5) 

^ G L2([0,oo),L2([0,l])), G L2([0, oo), iJi([0, 1])). 

Moreover, for any fixed T > 0, there exists a positive constant C (T) > which depends only on 
T and \\[uq,Vq)\\jj2 such that 

v^ix,t) >C{T). (3.6) 
Based on Theorem 1.3 (i) and Lemma 3.1, the following L^-convergence rates can be proved. 

Lemma 3.2 (L^-Convergence rates) Assume that the assumptions listed in Theorem 1.3 are 
satisfied. Then there exists a positive constant C , independent of e, such that 



f 

Jo 



dx 



dxdr < Ce^ 



(3.7) 



and 



dx 



dxdr < Ce. 



(3.8) 
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Proof. Set 



(3.9) 



Then wc deduce from (1.1), (1.5), (1.6) and (1.8)-(1.11) that {ijj^,6^) {x,t) satisfy the following 
initial-boundary value problem: 



with initial data 



and the boundary conditions 

(V^^e|)(o,^) = (v^e|)(l,^) = (o,o), t>o. 



(3.10) 



(3.11) 



(3.12) 



By the similar method to those in Subsection 2.2, we can obtain (3.7) and (3.8). We note that 
the convergence rate in (3.7) is the same as that shown in (2.48). However, the rate in (3.8) 
is slower than one obtained in (2.49), which is caused by the boundary effect. In fact, unlike 
Lemma 2.5, the regularity on solutions to the initial-boundary value problem (1.8)-(1.11) can 
not be improved. Consequently, the term J15 = — 2£ /q Jq tp^^u^^dxds (see (2.59) above) can 
only be controlled as follows: 



-2£ / / i^la^uLdxds < Ce [ f (ulSdxds + ej f {i^l^f dxds 
Jo Jo Jo Jo ^ ' Jo Jo 

<Ce + e / {iplJ^dxds. (3.13) 
Jo Jo 



Finally, based on Lemma 3.2, we can prove Theorem 1.3 (ii). In fact, using Sobolev inequality, 
we also have from (3.7) and (3.8) 



(u^ - u^) (t) 



L°°[0,1] 



< C 



(u'-u^) {t)f\\(u'-u°)jt) 



< C£4 



(3.14) 



and 



L°°[0,1] 



< C 



(v^ - v') it) 
This completes the proof of Theorem 1.3. 



(3.15) 



4 Appendix: Derivation of models 



In this section, we give the derives of conservation laws (1.1) and (1-3) for the convenience of 
the readers respectively, cf. [14, 25, 5]. 
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As in [14], one considers the case where /(c) = ac (a > 0) in the Keller-Scgcl model 
(1.20), which means that oxygen is consumed only when cells (bacteria) encounter the chemical 
(oxygen). The crucial step in the derivation is to make a change of variable through the Hopf- 
Cole transformation 

V = -C~'^Cx = -(Inc)a: (4.1) 

which was first introduced in [21] for a chemotaxis model proposed in [9] describing the chemo- 
tactic movement for non-diffusible chemicals (i.e. £ = 0), and was later applied in [12, 13] to 
study the nonlinear stability of traveling wave solutions. It turns out this transformation also 
extends its capacity to the full Keller-Segel model (1.20) for e / 0. Indeed, with the Hopf-Cole 
transformation (4.1), they derive the following viscous conservation laws from (1.20) 

ut + [eu^ -av)^ = euxx, 

(4.2) 

Substituting the scalings 

i=at, x = V = .f^v, -D = -, £ = - (4.3) 

into (4.2) and dropping the tildes for convenience, we obtain the system of conservation laws 
(1.1). By the transformation mentioned above, we get conservation laws (1.3) from (1.21). 
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